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J3 , Interferometric observations with three telescopes or more provide two 
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observables: closure phase information together with visibilities measure- 
ments. When using single-mode interferometers, both observables have to 
be redefined in the light of the coupling phenomenon between the incoming 
<^ ' wavefront and the fiber. We introduce in this paper the estimator of both 

in : 

. so-called modal visibility and modal closure phase. Then, we compute the 

^ . statistics of the two observables in presence of partial correction by Adaptive 

^ ■ Optics, paying attention on the correlation between the measurements. 

O ■ We find that the correlation coefficients are mostly zero and in any case 

i-^l never overtakes 1/2 for the visibilities, and 1/3 for the closure phases. From 

^ I this theoretical analysis, data reduction process using classical least square 

^ ! minimization is investigated. In the framework of the AMBER instrument, 

. the three beams recombiner of the VLTI, we simulate the observation of a 

single Gaussian source and we study the performances of the interferometer in 
^ ' terms of diameter measurements. We show that the observation is optimized, 

I i.e. that the Signal to Noise Ratio (SNR) of the diameter is maximal, when 

the full width half maximum (FWHM) of the source is roughly 1/2 of the 
mean resolution of the interferometer. We finally point out that in the case 
of an observation with 3 telescopes, neglecting the correlation between the 
measurements leads to overestimate the SNR by a factor of \/2. We infer that 
in any cases, this value is an upper limit, (c) 2008 Optical Society of America 

OCIS codes: 030.6600, 030.7060, 070.6020, 070.6110, 120.3180 



1. Introduction 

Thanks to the simultaneous recombination of the light arising from three telescopes, interfer- 
ometers such as I0NIC-3T on lOTA^ or AMBER on the VLTP are providing closure phase 
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measurements together with the modulus of the visibihty. Retrieval of phase information al- 
lows to scan the geometry of the source, hence opening the era of image reconstruction with 
infrared interferometric observations. However, with the current number of telescopes avail- 
able, direct image restoration requires many successive nights of observations.^ Thus, in most 
of the cases where the {u, v) coverage spans a relatively small number of spatial frequencies, 
the measurements have still to be analyzed in the light of model-fitting techniques. 

Furthermore, together with partial correction by Adaptive Optics (AO), many of the up- 
to-date interferometers are making use of waveguides that spatially filters the atmospheric 
corrugated wavefront, changing the turbulent phase fluctuations into random intensity vari- 
ations^ . The estimators that describes the visibility and the closure phase measurements 
obtained with such interferometers have to account for the coupling between the partially 
corrected wavefront and the fiber. From these appropriate estimators one can derive the 
statistical properties of the observables and properly investigate the performances of single- 
mode interferometers. 

In Section |21 we recall the spatial filtering properties of the waveguides in terms of inter- 
ferometric signal and we define the estimators of both the modal visibility and the closure 
phase. We investigate in Section El the covariance matrices of the observables with respect 
to atmospheric, photon and detector noises, paying particular attention to the correlation 
coefficients. Then, defining in Section 0] a general least square model fitting of the measure- 
ments, we analyze in Section El the ability of fiber optic interferometers to measure stellar 
diameters. 

2. Principles of fiber optic interferometry 

A full analysis of the signal arising from fiber optic interferometers has been theoretically 
described by Mege^ and summarized by TatuUi et al^ . We only recall here the important 
points for this paper, focusing on the coupling phenomenon between the incoming wavefront 
and the fiber, and on the observables that can be obtained from such interferometers. Figure 
H sketches the principle of a fiber optic interferometer, and reports the main technical terms 
that will be used all along this paper. 

2. A. Spatial filtering 

Introducing waveguides to carry/recombine the light in an interferometer allows to perform a 
spatial filtering of the incoming wavefront. It means that the phase corrugation of the wave- 
front are changed into intensity fiuctuations. In other words, the number of photometric and 
coherent (interferometric) photoevents at the output of the fibers depends on atmospheric 
fiuctuations and results in the coupling between the turbulent wavefront and the fibers^''' . 
Hence, spatial filtering can be seen as coupling coefficients, i.e. the fraction of (respectively 
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photometric and coherent) hght that is captured by the fibers. Such couphng coefficients are 
mathematically described by the following equations^'^ : 

p^(K) = Po(K*n/=o (1) 



P.,(K) = Po(K*r^);=;^^. (2) 

where is the visibility of the source and and T*-^ are resulting in respectively the auto- 
correlation and cross-correlation of the aberration-corrupted pupil weighted by the fiber 
single mode^'^° : 

~ /|P.(r)|Mr 
. ^ / P.{r)P;ir + Xm{r)^;ir + A/)dr 

J P,{r)P*{r)dr ^ ' 

where Pi{r) is the pupil function of the i^^ fiber optic telescope and ipi{r) is the aberration- 
corrupted wavefront incoming on the i^'^ pupil. and T*-' are respectively called the pho- 
tometric and interferometric peaks. The inverse Fourier transform of T* is called the photo- 
metric lobe (or antenna lobe as commonly named in radio) and the inverse Fourier transform 
of T^^ is called the interferometric lobe, po is the optimum coupling efficiency fixed by the 
fiber core design^^ . 

In the case where the source is unresolved by a single telescope (i.e. is much tighter than 
the photometric lobe), Eq. ^can be simplified: 

P^{V.)=Po J T{f)df = poS' (5) 

where S is the instantaneous Strehl ratio^^ . Moreover, if the visibility is constant over the 
range of the high frequency peak T*-', the interferometric coupling coefficient has also a 
simplified expression: 

\p^,m'=plS^S^\VM^,)\' (6) 

Under these conditions, the effect of spatial filtering by the fibers in the interferometric 
equation and, as a result, in the observables, is entirely characterized by the instantaneous 
Strehl ratio statistics. 

2.B. Estimation of the modal visibility 

We refer to TatuUi et al^ for a more complete description of the modal visibility. Note however 
that the expression of the coherent fiux at the spatial frequency fij is given by: 

/(/.,) = v/^p.,(K) (7) 
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where Ki and Kj are the number of photoevents of telescopes i and j before entering the 
fiber. An estimator of the modal visibility in the Fourier space is given by the ratio of the 
coherent flux by the photometric ones, assuming the latter are estimated independently 
through dedicated photometric outputs (see Fig. [T)): 



*\ kikj ^ 

r is the fraction of light selected for photometry at the output of the beam-splitter, and ki, 
kj are the photometric fluxes (after the fibers). 

2. C. Modal bispectrum and closure phase 

By definition, the closure phase is the phase of the so called image bispectrum B^im- The 
latter consists in the ensemble average of the triple product < I{fki)I{fim)I*{fkm) >■ It can 
be expressed from Eq. [7| as: 

Bklm =< KkKiKrnPkl{V,)pUy.)plm{y*) > (9) 

= KkKiKmpl jjj K(/)v;(/>;(/"). 

< T^\fki - /)r'"^(/M - f')T'"^'ifki - /") > d/d/'d/" (10) 

Using Roddier's formalism^'^ that demonstrated the bispectrum analysis to be a generaliza- 
tion to the optical of the well known phase closure method currently used in radio interfer- 
ometry, it is straightforward to notice that the quantity < T^''{f)T'''^'^^(f')T''"''* (f") > is non 
zero if f" = f + f and that in this case: 

< T'\f)T'^[f')T'^\f") >= (11) 



KkKiKMf)VM')V:U-") = BM\ /') (12) 

where iV(/, /') is proportional to the overlap area of three pupil images shifted apart by the 
spacings /, /' and / " = / + /', and B^(f,f') is the bispectrum of the source. Hence the 
modal bispectrum can be rewritten: 

B,,^ = pljj BM, /) ^^^'''^{(^p^'^ d/d/ (13) 

Thus, the modal bispectrum arising from fiber optic interferometers is the source bispectrum 
integrated over the overlap area N{fki — f,fim — f')- Hence, as the modal visibility does 
not equal in general the object visibility, the modal bispectrum does not coincide with the 
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source bispectrum. Nevertheless, if the source spectrum is constant over the overlap area 
N{fki — /, fim. — f')i Eq. Uni takes a simphfied form: 

B,im. = plB.ifku fim) II ^^?[^d/d/' (14) 
In this case, the modal bispectrum is proportional to that of the source. 
3. The covariance matrices 

We propose to characterize the statistics of the square visibility and that of the closure phase 
(i.e. the bispectrum phase) by computing their respective covariance matrix. Our objective 
is twofold: derive the error associated to each observable and investigate the degree of de- 
pendency of each observable through their correlation coefficients. In order to do so, we use 
the spatially continuous model of photodetection introduced by Goodman^^'^^ where the 
signal is corrupted by three different types of noise: (i) the signal photon noise; (ii) the ad- 
ditive Gaussian noise of global variance which arises from the detector and from thermal 
emission; (iii) the atmospheric noise resulting from the coupling efficiency variations due to 
the turbulence. To simplify the calculations, we assume that the source is unresolved by a 
single aperture, such that the low frequency coupling coefficients verify Eq. El and that the 
source visibility is constant over the range of the interferometric peaks, such that the high 
frequency coupling coefficients verify Eq. IHl These assumptions drive to neglect the modal 
speckle noise regime - it has been shown that the modal speckle noise is rejected towards 
negative magnitudes and only affects very bright sources® - and to only focus on "photon 
noise" and "detector noise" regimes . The full calculations of the covariance coefficients of 
the visibilities and of the closure phases are done in Appendix A. They lead to relatively 
complicated formulae which depend on the Strehl statistics. Using a simple analytical ap- 
proach, we derive in Appendix Bthe mean and the variance of the Strehl as a function of 
the turbulence strength and the level of AO correction. The relative error of the Strehl is 
bounded between two limit values^^'^^ : 

(Perfect correction) < — < 1 (No correction) (15) 

S 

Table ^ gives the expressions of the limiting values of the variance of the visibility and the 
closure phase for a point source in both "photon noise" and "detector noise" regimes. The 
error of the visibility will be deeply used in the next section to derive the performances of the 
Very Large Telescope Interferometer (VLTI) with regards to single sources diameter measure- 
ments. Let us concentrate in this part on the correlation coefficients of the visibilities and 
the closure phases, respectively. Their limiting values are summarized in Table |21 Clearly, for 
visibilities the correlation coefficients are null in the photon noise regime when no telescope 
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is in common (see Fig. |21) and are always smaller than 1/2 otherwise, and for the closure 
phases they are null when no baseline is in common (see Fig.lHl) and are always smaller than 
1/3 otherwise. Furthermore, the relative number of null elements in the covariance matrices 
rapidly increases with the number of telescopes. It implies that if there are (A^ — 1) (iV — 2)/2 
linearly independent closure phase relations"'^^ , the whole set of A^(A^ — 1)(A^ — 2)/6 closure 
phase relations can be considered in first approximation as statistically independent. 

It would be of great interest to pursue this study by comparing the statistics of the 
closure phase to that of the phase measured by phase referencing technique. The formalism 
presented in Appendix A can indeed be transposed to the phase referencing case, but it 
requires a dedicated analysis which is beyond the scope of the present paper. 

4. Model fitting 

The present generation of interferometers only provides small number of telescopes (basically 
2 or 3). In such a case, it is frequent that the lack of spatial frequencies in the {u, v) coverage 
prevents from image reconstruction of the studied object. Hence, the measurements have to 
be analyzed in the light of a model of the object. We propose here a simple model-fitting 
of the observables. 

Let us define Oe as the normalized spectrum of the object model characterized by a set 
of parameters 6. From this spectrum we derive the model of the observables that have to be 
fitted according to the measurements: 



Square visibility 
Closure Phase 



\0s 



= arctan 



lm{Ge 
Re{Gg 



where G^'"^ = Og{fki)Oe{fim)Og{fkm) is the bispectrum model. Then, the estimated pa- 
rameters 6 constraining at best the observations are obtained by minimizing the distance 
between the model and the measurements. Assuming Gaussian statistics for the observables, 
the distance is given by the well known x^'- 



2 _ 2 I 2 
Xtot — Xv^ + Xcf,], 



with 



- \0g 



^-1 



V^-\Oe\ 



(16) 
(17) 



X(j)B 



where X'^ denotes the transpose of the vector X. Equation supposes that the two ob- 
servables of different nature (i.e. the visibility and the closure phase) are uncorrelated. Such 
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assumption seems quite reasonable since we have shown that two measurements of same 
nature are already poorly correlated. 

It can be argued that, in order to avoid phase discontinuities issues, it might be better to 
work on the phasor (i.e. the average bispectrum) than on the phase of the average bispectrum. 
In that case the constraint takes the form: 



Im BGl 



(7-1 

lm{BG;] 



Im BGl 



(19) 



Such constraint does not appear to be appropriate for two reasons: (i) constraining the 
average bispectrum or the closure phase of the average bispectrum is strictly equivalent 
when the latter shows good SNR; (ii) the covariance depends on the model which makes 
the minimization subject to biases due to improper noise estimates. For very noisy data, 
the closure phase shows a lot of discontinuities and its probability law, wrapped around 
[— 7r,7r], tends towards a uniform law. Such a case, for which the fitting method is no 
more optimal, neither for the closure phase nor for the bispectrum, corresponds in practice 
to the sensitivity limit of the instrument. 
Putting Equation ^1 into the generic form: 

x'{e) = MoCj^Mj (20) 

the error of the estimated parameters 6 writes^'^ : 

a\9) = Dmg{[A^C:^A]-'} (21) 



where 

A- 



(22) 



89 

5. Applications 

5. A. Observing a single Gaussian source 

We propose in this section to simulate the observation of a single Gaussian source of FWHM 
(To, with AMBER, the three beam recombiner of the VLTl^ . Since the object is centro- 
symmetrical, closure phase is not relevant and hence all the information is contained in the 
visibility alone. For sake of simplicity, we first neglect the contribution of the correlation 
coefficients. Their effect will be studied in Section 15. CI We adopt the standard instrumental 
configuration of AMBER^^ in which the signal is sampled over Npi^ = 16 pixels. We choose 
a spectral resolution of 35 in the K band (2.2/im), an integration time of 30ms per inter- 
ferogram, a transmission coefficient r = 0.5, a readout noise of 15e~/pix and an optimized 
coupling coefficient of po = 0.8^^ . We observe an object with 2 Unit Telescopes (UT2 and 
UT4, D = 8m) during 4H (half of the time on the object, half of the time on the calibrator) 
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with 5min per frequency point which, together with the integration time per interferogram, 
leads to 5000 samples per frequency point. Note that all along the observation, the length of 
the projected baseline is quite constant with B = 45m ± 2m. Finally, we assume a turbulence 
strength of D/vq = 5 and a typical AO correction oi S = 0.5. 

Fig. m shows the SNR of the object size as the function of the magnitude. We consider 
different sizes: (Tq = Imas, (Tq = 3.8mas, cTo = 5.4mas, and (Tq = lO.Smas, which are to 
be compared to \/B ~ lOmas, the resolution of the interferometer. As expected, a general 
profile can be seen with two well known regimes: the "photon noise" regime for bright 
sources and the "detector noise" regime for faint sources. Defining the limiting magnitude 
as the magnitude for which the SNR is equal to 1, we find /C ~ 11 — 13 according to the size 
of the source. 

Clearly, the SNR first increases and then decreases with the size of the source, reaching a 
maximum around (To = 4mas. This phenomenon can be understood as follows: for marginally 
resolved sources, the parameters of the fit are barely constrained and the SNR is small. It 
increases with the size up to a point where the available projected baseline range does not 
match anymore the frequency content of the object. From there, the SNR begins to drop. 
This trend stands for all observing conditions, but as we show in the next section the exact 
position of the maximum depends on the noise regime. 



5.B. Optimizing the baseline 

We simulate different configurations with two telescopes that span the range of average 
projected baseline, respectively: (a) B = 45m (+); (b) B = 56m (x); (c) B = 83m (o); (d) 
B = 100m (A); (e) B = 124m (□). The declination of the source is arbitrarily set to —25°. 
The source is supposed to be observed between — 3h to 3h from the zenith. The parameters 
of the different configurations used are summarized in TableEland corresponding (m, v) plane 
coverages are shown in Fig. El (up). 

We compute the SNR of the diameter for all these configurations, in both "detector noise" 
{K = 13) and "photon noise" {K = 2) regimes, respectively. In the "detector noise" regime, 
the error of the visibility is independent on the diameter, and since in this case Eq. |^ 



tells us that a (6) is directly inversely proportional to 



dM, 



dd 



the SNR of the diameter is 



maximum when the product of the diameter by the derivative of the model is maximum too, 
namely when: 
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This leads to: 



7r 



V2 <B> 



0.38 



A 



<B> 



(23) 



(24) 
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where < 5 > is th average projected baseline. 

In the "photon noise" regime, the error of the visibihty is strongly dependent on the size. 
In that case, the maximum of the SNR occurs when the projected baseline range scanned 
by the interferometer matches the object frequency content. Considering the frequency for 
which the visibility is 1/e, it comes: 

2Jln(2) A A 
ao = ^ ^ ' ^ - 0.54—— 25 

Tx < B > < B > 

Fig. El (middle and bottom) illustrates these results. The SNR of the diameter is plotted as 
a function of the size in \/B units, for the five geometrical configurations selected above. 
We can see that the whole curves present the same behavior, and especially the same max- 
imum. This maximum verifies Eq.'s and 1^ for the "detector noise" and the "photon 
noise" regimes, respectively. We conclude that the mean projected baseline optimizing the 
observation of an object of diameter o"o is given by: 

< S >^ 0.5— (26) 

Furthermore we stress that the performances degrades rapidly when this criterion is not 
respected. 

5. C. Effect of the correlation coefficients on the error bars 

We have shown in Section El that the correlations between the visibilities could reach a 
maximum value of p = 0.5. Since the estimator arising from maximum likelihood (i.e. the 
minimization) is unbiased, the expected value of the fitted parameters does not depend on 
whether the correlation coefficients are introduced or not. We analyze here the effect of the 
correlation coefficients on the derived error. 

We simulate the observation of a Gaussian source with three UTs, together with the 
instrumental parameters of Section El and we choose cTq = 0.5A/ < B > for the FWHM of 
the source. Fig. El shows, in both pure turbulent (cr^ = S) and fully AO corrected (as = 0) 
cases, the SNR of Uo, with and without taking into account the correlation coefficients. 
It results that considering uncorrelated measurements leads to underestimate the error (or 
overestimate the SNR) by a factor of a/2, roughly. Given that the ratio of the null elements 
versus the non null elements in the covariance matrix increases with the number of telescopes, 
we infer that this factor -\/2 is an upper limit. 

6. Summary 

We have computed the theoretical covariance matrices of the modal visibility and the modal 
closure phase in the presence of partial AO correction, when the measurements are cor- 
rupted by atmospheric, photon and detector noises. In the photon noise regime and for an 
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interferometer with a large number of telescopes, the measurements are most of the time 
uncorrelated. In any case, the correlation coefficients are always smaller than 1/2 for the 
visibilities and smaller than 1/3 for the closure phases. 

Then from a classical least square approach, we have investigated the ability of interferom- 
eters to measure stellar diameters. In the light of the AMBER experiment we have found a 
limiting magnitude in the range K = 11 — 13 depending on the size of the source. At last, we 
have shown that the observation is optimized when the mean resolution of the interferometer 
is equal to twice the stellar diameter. 

APPENDIX A: STATISTICS OF THE OBSERVABLES 

1. General formalism 

In order to compute the moments of the spectral density, we use the spatially continuous 
model of photodetection process of Goodman^^ where the detected signal is corrupted by 
photon noise, by additive Gaussian noise e of variance and by the turbulent atmosphere. 
It takes the form: 

K 

s{x,y) = ^5{x - Xn) +e{x) (A-1) 

n=l 

and its Fourier transform: 

K 

S{f) = exp (-2m/x„) + e(/) (A-2) 

n=l 

2. The modal visibility 

Such calculation has been already done in TatuUi et al^ . We give here the results assuming 
that the low and high coupling coefficients verify Eq. |31 and El respectively. We remind that 
following equations assume also that the telescope transmissions are all equal, i.e. Ki = 
K/Nteh where Ki is the number of photoevents coming from telescope i and Ntei is the 
number of the telescope, and that the level of corrections of the Adaptive Optics systems 
are the same for all the telescopes. 

The square relative error of the modal visibility can be seen as the sum of two contributions: 

= £l{K, S) + £l{K, a\ S) (A-3) 

with the photon noise square relative error and E\ the additive noise relative error, as 
described in Tatulli et al^ . The same way, the covariance can be cut in two terms: 

""^1:;; = Cp(K, S) + Ca{K, cy\ S) (A-4) 
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Here we must pay attention that two cases can occur (see Fig. |21 in the body of the text): 
(i) either the basehnes fij and fki come from two distinct pairs of telescopes (ii) either the 
basehnes have one telescope in common, let say j = k, which drives to extra correlation 
between the visibilities. 

The results, i.e. the diagonal and non diagonal terms of the visibility covariance matrix, 
are summarized in table El 

3. The closure phase 

The estimator of the bispectrum between telescopes k is defined by: 

Qi,k = m,)Sifjk)S*{U) (A-5) 

Chelli^^ has shown that the covariance (per sample) on the closure phase does only depend 
on the modulus of the bispectrum, and hence assuming a centro-symmetrical source, could 
be written: 

Gov = l^^^j^j^^f^ (A-6) 

with 

E(Q,,fc) = <Z'?(/,,)T(/,fe)?(/fcO (A-7) 

(A-8) 



where i{f) is the normalized spectral density. 

For sake of simplicity we analyze the second order moments of the estimator in two cases 
separately: the "photon noise" case and the "detector noise" case. Then the statistics with 
regards to the atmosphere are taken into account. 

4- The photon noise case 

Following Goodman's formalism^^'^^ the second order moments are: 



^K{\Qijk\^} 



+K 
+1^ 



(A-9) 



K 



+ 2\i{uMUok)? + 2\i{uM{Uk)? + 2|?(/,,,)rR/.fc)p 



+2K' 

+ K 



i{fij)i{fjk)i*{fik) 



(A-10) 
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For the covariance, two cases may appear (see Fig. El in the body of the text): (i) both 
triplets of telescopes k and l,m,n) are different or (ii) two telescopes are part of both 
triplets, in that case there is a baseline in common, let say fij = fi^. This add extra terms 
in the covariance. 



EKiQijkQlmn] = K 



'iUij)'i'UikyUik)i\flm)i*Urnn)iUln] 



i-f fij flm 



+K 



iifknWikKifln) + iifkn)iifjk)i*if 
iifjk)i*ifik)i*{fmn)i{fln) 



+K''\i{fij)\^ i{fkn)i{fik)i*{fin) + i{fkn)i{fjk)i*{f 



ExiQijkQlmn} = K i{fij)i{fjk)i*{fik)i{flm)i{fm.n)i*{fln) 



mn / 



(A-11) 



if fij = fi 



ij J Im 

'mmn)\'\Uk)\' + mmn)\'mk)\' 

iifjkKifikKifmnWln) 
^{fiJWJkr{f^kmfm.n)\' + IKflnW) 
KflmWmn)i*{fln){\i{fjk)\^ + R/.fc)!' 



(A-12) 



5. The detector noise case 

We suppose a zero mean Gaussian detector noise of variance a^. 



MUjk\'} = K m,)?\i{fjk)\Mf 



lik 



^if^^)mfJk)\' + m,)\'mk)\' + m,k)\'mk)\' 

2 



lik 



(A-13) 
(A-14) 



^e{QijkQ*lmn] = ^ \i{fij)i{fjk)i*{fik)i*{flm)i*{fmn)i{fln) 
if fij flm 



-Na^K 



Kfij)i*{fik)i*{fmn)i(,fln) 



2^6 



(A-15) 
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'E^eiQijkQlmn} = K i{fij)i{f jk)i\fik)iUlrn)i{frnnY [fin] 



IhNa 
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(A-16) 



Once taking into account the statistics with regards to the atmospheric turbulence, the 
final expression are given in both "photon noise" and "detector noise" regime in tables El 
and El respectively. 

The expression of the covariance in the general case is obtained by adding the covariance 
in each regime: 



Gov ^Imn} = CoVphot {^ijk, ^Imn} + Cov^et {^ijk, ^Imn} 



(A-17) 



APPENDIX B: MODELLING PARTIAL ADAPTIVE OPTICS CORRECTION 



We consider a simplified model where the long exposure AO corrected transfer function can 
be divided in two weighted components: one perfect transfer function and one turbulent 
transfer function where the weight h describes the strength of the correction: 



< r (/) > 



[h + {i-h)l3^if)], he [0,1] 
1 



exp 



(B-i: 



where T>^{f) is the AO corrected structure function. For sake of simplicity we assume that 
the transfer function of the turbulence B^{f) is Gaussian, hence: 



= exp ( ) ,as 



2 rp 
6.88 A 



(B-2) 



Going further we can notice that the AO corrected transfer function of the atmosphere can 
be decomposed in a low frequency term (halo) and a high frequency term (coherent energy) 



exp 



-2^<a(/) 



and that, 



Hence we have: 



exp 



h + {l-h)B^{f) 



(B-3) 
(B-4) 



h = exp(-aj) = E,(Z}/ro, N,) (B-5) 

where the coherent energy E^ depends on the level of correction (number of Zernike A^^) and 
the strength of the turbulence D/tq. 
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Using the definition of tlie long exposure Strelil ratio: 

S = J <r(/)>d/ (B-6) 

it comes: 



S = h + {1- h)-^f-^ = E, + ^ (B-7) 



JO 

The second order moment of the instantaneous Strehl ratio writes: 



= II <Tif)T{f)>dfdf 

Pi{r)Pi{r + p)Pi{r)P,{r + p) 
X < ^'i(r)^'*(r + p)-^*{r')'^i{r' + p)> drdpdr'dp' (B-8) 

where Pi{r) is the pupil function and \l/j(r) is the complex amplitude of the AO corrected 
wavefront. We use Korff's^^ derivation of the moments of the complex amplitude of the 
wavefronts to introduce in the equations linear combinations of the structure function {T>{r)) 
at different spatial frequencies. Assuming the structure function to be stationary, it comes: 



52 = P,{r)P,{r + p)P,{r')P,{r' + p). 

exp {-\[V{p) + V{p) + V{r' - r) + V{r' + p - r - p)]] 
exp { — |[r'(r' + p' — r) + V{r' — t — p)] } 



drdpdr dp (B-9) 



Using Eq. IB-41 previous expression can be analytically developed for /i = E^ =]0.5, 1], i.e. for 
good AO correction levels^ . However we know that, without AO correction and assuming 
circular Gaussian statistics for the complex amplitude of the wavefront, we have as = S^^'^"^ 
. Hence we can perform an interpolation between the good AO corrections and the pure 
turbulent cases. Fig [7| shows the SNR of the Strehl ratio as a function respectively of the 
number of Zernike (i.e. the correction level) and the average Strehl ratio S, at different D/tq. 
For each curve, we point out where from the interpolation has been done. 
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Table 1. Variance of the visibility and the closure phase for a point source. Npi^ 
corresponds to the number of pixels that sample the interferogram. Npi^ > 2 
is required to respect the Shannon criterion, a is the detector noise per pixel. 



Observables 



Variance of the observables (point source) 



Photon noise regime {K ^ 1) 



Full AO correction 



No AO correction 



Detector noise regime 
(Z< 1) 



y2 



(l~r) 



PqSK 



(1-r) + r 



PoSK 



(l-r)4 



Nf , 

tel 



PoSK 



PoSK 



3^ 



P^S K 
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Table 2. Correlation coefficients of the visibilities and the closure phases. For 
the visibility, two cases are considered: one telescope is common to both base- 
lines, hence the triplet of telescopes is forming a triangle (a so-called closure), 
or not (see Fig. |21). For the closure phases, two cases are investigated as well: 
one baseline belongs to both closure phases or not (see Fig. E}. 







Correlation coefficient p (point source) 




Observables 


Photon noise regime {K 3> 1) 












Detector noise regime 










(Z< 1) 






Full AO correction 


No AO correction 






if closure 


i 

< - 
2 


i 

< - 
4 






otherwise 








3 

3 + Npix 




if baseline in common 


1 
3 


1 

< - 
6 


3 














otherwise 
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Table 3. Description of the considered telescopes configurations. Are given 
the name of the telescopes pairs, their respective average projected baseline as 
well as the related symbols used in Fig. The telescopes are the four 8 meters 
Unit telescopes (UTs) of the VLTI. The declination of the source is arbitrarily 
set to —25°. The source is supposed to be observed between — 3h to 3h from 
the zenith. 



Symbol 


Telescopes 


Average projected baseline 


+ 


UT2-UT3 


45m 


X 


UT1-UT2 


56m 


o 


UT2-UT4 


83m 


A 


UT1-UT3 


100m 


□ 


UT1-UT4 


124m 



19 



Table 4. Description of the elements of the visibility covariance matrix. 



Obscrvables 



Covariance cocfBcicnts 



Visibility 



Diagonal 



coefficients 



(1- 



r)V?(/,i)5^ 



(l-r)2v;2(/,j)5^ 



2(4 +-g) 



Ntel 

PoK 



(1 
2Na^ 



(i-r)2v;2(/,,)s' 



(i-r)4y4{/,,)s* 



t25 



pIk' (1 



Non diagonal 



coefficients 



(1- 

if j: 



(l-r)2K2(/,,)y,2(A,)5" 



1 Ntel 
tS poK 



{I - r)v,Hfij)v.Hfki)s 



Nil 

pW 



pW' 



PoK 



vKUl) 



(l^r)-^V^{fij)Vnfkl)S^ 



Nil 

pW' 



Nil 

pW' 



Table 5. Description of the elements of the closure phase covariance matrix, in 
the "photon noise" regime. 



Obscrvablcs 



Covariance coefficients 



Closure phase 



Cov{$ijfe,'I>;„„} 



Diagonal 



coefficients 



Non diagonal 
coefficients 



mUij)?\vAfjk)?\v.Ui}^)\'s'' 



m{.hjmv.{fjk)\^\vAfik)\^s'' 



pW' 



- 245] 



2!v.(/i,)|2|v;(/,.fc)|2|K(/.fc)P5;^ 



pW 



Ntel 

PoK 



Ntel 

PoK 



if .Uj + fl, 



|K(/fcn)l'[A^tei5 +245] 



2V^ifij)V*ifjk)\V:(fik)V:iflm)V:ifmn)\V^Ul^)S° 

[\V*{f^j)\^\VAfkn)\^]S* + [V;(/fcn)V7(/mn)V7(/^fc) + K (/fcn ) (/ifc ) (/in )] [A^M^" + 34^^ 



Nil 



'2V^(hj)VMjk)\V:ihk)V:iflrr,)V:ifmn}\V^ifln)S° 



pW 



[V*{fjt:)V:{f,k)V:{fj„)V,{fin)][NtelS^ +4alS^] + [14 (/fcn) (/„) (/ifc ) + {/fc„ )K (/j fc ) (/jn )] I ) | ^ [^^^ + 45"] 



2 y, if, J )V*{fjk)\ V* Uik)V* Ulm)V* {fmn)\V4flr,)S'' 

[v.{fij)VAfjk)v:(fik)m{fi,.)\^ + |K(/i„)p) + v*(/,,)V4/in)K*(/in)(|v;{/,fe)P + \VAfik)ms' 



-air 



2v;(/y)K(/,fc)|y^*(Afe)K*(/im)v;*(/,„„)|K(/,„)5° 



Ntel J. 

PoK 



Ntel 

PoK 

flm 



Table 6. Description of the elements of the closure phase covariance matrix, in 
the "detector noise" regime. 



Obscrvablcs 



Covariance coefficients 



Closure phase 



Gov {'!>ijk,'S>lmn} 



Diagonal 



coefficients 



Non diagonal 
coefficients 



2|V;{/,,)P|K(/,fe)P|K(/ife)P5° 



mUij)?\v.(fjt.)?\vAfi>.W'^ 



if f^J 7^ fl, 



[v*{fjk) v: ihk) v: {fj^)v. (fir,)] Na^s 



2K )v.ifjk)\v:(hk) v: UimW: {fm„) 1 )5° 



pW 



if fij — flm 
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Photometric output 




Interferometric output 



Pliotometric output 



Fig. 1. Sketch of a fiber optic interferometer. Thanks to single mode fibers, 
beams coming for each telescope are spatially filtered. Then for each beam, 
a beam splitter is selecting one fraction of the light in order to estimate the 
photometric fiux through dedicated photometric output. The remaining part 
of the light is recombined with the other beams (recombination is seen here 
as a black box, and the recombination mode - spatial or temporal - is not 
a relevant parameter in this study). At the output of the recombination, the 
coherent flux is estimated thanks to the so-called interferometric output. 
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Fig. 2. This figure shows the configurations that have to be considered when 
computing the covariance of the visibilities. Respectively no telescope in com- 
mon (left) or one telescope in common, hence forming a closure (right). 
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Fig. 3. This figure shows configurations that have to be considered when com- 
puting the covariance of the closure phases. Respectively no baseline (left) or 
one baseline (right) in common between both triplets of telescopes. 
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K — magnitude 



Fig. 4. SNR of the diameter cTo as a function of the magnitude. Different 
sizes are considered: o"o = Imas (dash-dotted hne) , cTq = 3.8mas (dotted 
hne), (Jo = 5.4mas (sohd hne), ao = lO.Smas (dashed hne). The resolution is 
\/B = lOmas. S = 0.5, D/ro = 5. 
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Fig. 5. Up: {u, v) coverage for different UTs configurations. The average pro- 
jected baseline length is respectively: B = 45m (+), B = 56m (x), B = 83m 
(o), B = 100m (A), B = 124m (□). More informations about the different 
configurations considered are given in Table El Middle: SNR of the diameter 
as a function of the size given in fraction of the interferometer resolution X/B, 
in the "detector noise" regime. Vertical line shows the size where the SNR is 

maximum. Bottom: same as above, but in the "photon noise" regime. 
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Fig. 6. SNR of the diameter as a function of the magnitude, for perfect AO 
correction (upper curves, as = 0) and in the pure turbulent case (lower curves, 
ag = S). SNR is computed assuming diagonal covariance of the measurements 
(no correlation: solid line) or considering correlation coefficients as computed 
in Section ini (dashed lines). 
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Fig. 7. SNR of the Strehl ratio as a function of the level of correction in number 
of corrected Zernike (top), and as a function of the long exposure Strehl ratio 
(bottom). Results are given for D/vq = 1 (solid line), D/r^ = 2 (dashed line), 
D/tq = 5 (dash-dotted line), D/tq = 10 (dotted line). Vertical lines show from 
where the interpolation has been performed. 
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